Prom Slater to Mott-Heisenberg physics: The antiferromagnetic phase of the Hubbard 

model 
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We study the optical conductivity of the one-band Hubbard model in the Neel state at half filling 
at T = using the dynamical mean-field theory. For small values of the Coulomb parameter clear 
signatures of a Slater insulator expected from a weak-coupling theory are found, while the strongly 
correlated system can be well described in terms of a Mott-Heisenberg picture. However, in contrast 
to the paramagnet, we do not find any evidence for a transition between these two limiting cases 
but rather a smooth crossover as a function of the Coulomb interaction. 



I. INTRODUCTION 

The microscopic description of magnetism and metal- 
insulator transitions constitutes one of the major re- 
search activities in modern solid state theory. For ex- 
ample, transition metal compounds like V2O3, LaTiC-3, 
NiS2-xSe x or the cuprates show metal-insulator tran- 
sitions and magnetic order depending on composition, 
pressure or other control parameters^ One interesting 

and controversial question concerns the description of 
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the optical properties of these materials, 6 in partic- 
ular whether the fundamental physics is governed by the 
broken translational symmetry e.g. in the Neel state or 
rather by correlations^^ i.e. the formation of so-called 
Hubbard bands with an energy gap of the order of the 
relevant Coulomb repulsion. 

The simplest model showing both magnetism and a 
correlation induced metal-insulator transition (MIT) is 
the one-band Hubbard model^ 
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where c^J annihilates (creates) an electron at site i with 
spin <t, rii a = c\ a c i(J , tij denotes the hopping ampli- 
tude between sites i and j and U is the local Coulomb 
repulsion. Usually, one ignores longer range hopping 
processes and concentrates on nearest neighbor hopping 
only. Considerable progress in understanding the physics 
of this simple but nevertheless non-trivial model has been 
achieved in the last decade through the development of 
the dynamical mean-field theory (DMFT)^^. In par- 
ticular, the phase diagram for the Hubbard model on a 
simple cubic lattice with nearest-neighbor hopping is very 
well understood&Siifi. The major results are compiled in 
the schematic phase diagram in Fig.^ At half filling the 
physics is dominated by an antiferromagnetic insulating 
phase (AFI) for all U > 0. For finite doping, the antifer- 
romagnetic phase persists up to a critical doping <5,~2iAi 
and in addition shows phase separationiii^. For very 
large values of U the antiferromagnetic phase is replaced 
by a small region of Nagaoka type ferromagnetism^*i4 
An appealing property of the DMFT is the possibility 
to calculate transport quantities in a very simple fashion. 



FIG. 1: Schematic DMFT phase diagram of the Hubbard 
model with nearest neighbor hopping on a simple cubic lat- 
tice. PM denotes the paramagnetic metal, AFI the antiferro- 
magnetic insulator, FM the ferromagnetic metal and AF/PS 
phase separated antiferromagnetism. 



Due to the local nature of the theory, vertex corrections 
to the leading particle-hole bubble of the current-current 
correlation function vanish identically^S*^ i.e. one needs 
to calculate the bare bubble only. This has been exten- 
sively used to study the optical conductivity and various 
other transport properties in the paramagnetic phase of 
the Hubbard modelASiiii On the other hand, up to now 
a comparable investigation of the optical properties of 
symmetry broken phases, in particular the Neel state at 
half filling, has not been performed. However, such an 
investigation is interesting for several reasons. First, the 
insulating phase in real materials is in many cases accom- 
panied by magnetic or orbital ordering, typically of the 
Neel type. To what extent the model can describe 
the optical properties of ordered insulating phases has 



2 



up to now not been studied in detail. Second, it is well- 
known that the restriction of the Hubbard model to the 
paramagnetic state at half filling shows a metal-insulator 
transition2i±2iiLi& at a finite critical Z7 C > which is of 
first order iSiii It might be argued that for the Neel state 
a similar situation can occur. At small U a weak coupling 
theory is expected to give accurate results, leading to a 
band or Slater insulator 2 due to the doubled unit cell in 
the Neel state. At large U, on the other hand, the Hub- 
bard model is known to reduce to an effective Heisenberg 
modeiiS with localized moments from the onset. It is an 
open question whether these two limits are linked conti- 
nously or via a phase transition at some finite value of 
the Coulomb interaction U. 

The paper is organized as follows: In the next section 
the derivation of an expression for the optical conductiv- 
ity in the Neel state is presented. The results obtained 
for the optical conductivity of the Hubbard model with 
nearest neighbor hopping on a simple cubic lattice at half 
filling are presented and discussed in section III. A con- 
clusion and outlook finish the paper in section IV. 



II. OPTICAL CONDUCTIVITY FOR THE NEEL 
STATE IN DMFT 



In the Neel state, the DMFT equations have to be 
modified to account for two inequivalent sublattices A 
and B (see Fig. |3 left panel) with self-energies E" 4 ^ 
£B | 9 i 20 rp Q j.jjjg enc ^ we m t r0( j uce operators a\J and b\J 



where we introduced the spinors 
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FIG. 2: Left: Schematic view of the AB sublattice decom- 
position suitable for the treatment of the Neel state. Right: 
Magnetic Brillouin zone (MBZ, 1. BZ of the Neel state). 

which act on sublattice A and B respectively. In the case 
of nearest-neighbor hopping only, the kinetic part of the 
Hamiltonian 0} can then be written as 
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After Fourier transforming this expression we obtain 
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and ek is the dispersion on the bipartite lattice. The 
prime on the sum indicates that the summation is over 
all values of k in the magnetic Brillouin zone (MBZ) (see 
Fig. |21 right panel). Within this notation, the Green 
function becomes a matrix in the two sublattices, 



Scj 

-£k 



-£k 



(2) 



where C,^ B = z + fi — T,^ B ■ From now on we employ the 
symmetry = ^§ = ( a of the Neel state and drop the 
indices Aj B. Using this formalism, the current operator 
is given by 
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with Vk = Vkfik as usual. If we consider a lattice for 
which the conductivity tensor is diagonal, the elements 
an = a can be calculated from (D is the spatial dimen- 
sion of the lattice) 



D ■ cr(w) 
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with the current-current correlation function 
{{3l\3l))iv = e 2 ££'«, 

aa' kk' 

Again, due to the symmetry of the lattice, the index I can 
be dropped. The most important simplification arises 
from the locality of two-particle self-energies within the 
DMFTi 15 ! 20 ! 21 Note that in the present formulation the 
proper locality of the two-particle self-energies is still en- 
sured, because in the DMFT as defined by equ. (0 no 
dynamical correlations between the A and B sublattices 
are introduced. In analogy to the paramagnetic case this 
allows us to carry out the k sums in diagrams containing 
two-particle self-energy insertions independently at each 
vertex. Since the single particle propagators only depend 
on k through the even function ek and the are of odd 
parity, the sum over their product vanishes. As a result, 
we obtain the exact expression for the current-current 
correlation function in the DMFT, 

ul n a k 

x [(( a k CT ; a L»i"»+^ (( 6 k CT ; b D)^ n 

+ ((a ka ,blj luJn+w ((a k(j ; blj iUn 
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In terms of the Green function matrix elements in (j2J we 
can rewrite this as 
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Next, we convert the k sum into an energy integral by 
introducing the average squared velocity 
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III. RESULTS 

In the following we present results for the optical prop- 
erties of the Hubbard model on a simple hypercubic lat- 
tice with nearest neighbor hopping at half filling and 
T = in the DM FT. The hopping matrix element is 
chosen as t — t* /v4D, which ensures the correct scal- 
ing of the kinetic energy in the limit D — > ooi£ As en- 
ergy unit it is convenient to use the bandwidth W of the 
system at U = 0. Since the Gaussian density of states 
(DOS) of the simple hypercubic lattice in the limit of 
infinite spatial dimensions has no real band edges, we 
take W — At* as a reasonable value. Note that for 
this choice of W, (i) the spectral weight of the Gaus- 
sian is exhausted by 99% between ui = ±W/2 and (ii) 
the paramagnetic metal-insulator transition will occur at 
U c w At* = The effective quantum impurity model 

of the DMFT is solved using Wilson's Numerical Rcnor- 
malization Group method (NRG)^ suitably extended 
for dynamical quantities and spin-polarizationi^S The 
calculations were performed with a discretization param- 
eter A = 2, keeping 800 states. Dynamical quantities 
were calculated with a Gaussian logarithmic broadening 
of 0.6. Occasional checks with 1600 states or smaller A 
showed sufficient robustness of the results. 



Making furthermore use of the spectral representation of 
the Green functions, the frequency sum can be evaluated 
in a straightforward way and finally we obtain for the 
conductivity 



A. Single-particle properties 
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x [A a {e,uj')A s {e,uj' + uj) + B a {e,uj')B a {€,uj' + u:)] (4) 



with 



and 



7T 



B a (e,u) = <SmG£ B (e,uj + i5) . 

Here /(w) is the Fermi function and c collects various 
constants. Note that the form (0} is reminiscent of the 
results found in the case of superconductivity, which is 
discussed at length e.g. in the book by Mahan.2 2 . Conse- 
quently, one can expect to obtain similar features from 
the evaluation of QJ. 

In order to proceed with the calculation, it is necessary 
to specify the actual lattice structure and the correspond- 
ing non-interacting dispersion in equation ©. For the 
hypercubic lattice^ (f 2 )e oc p^(e) is a simple Gaussian, 
and the integration over e can be performed analytically. 
For the details of this calculation see Appendix A. 




FIG. 3: Spin-resolved DOS for the antiferromagnetically or- 
dered phase of the half-filled Hubbard model for different val- 
ues of U/W. The inset shows the magnetization as function 
of U/W. 

Before discussing the optical conductivity calculated 
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from it is helpful to review the single particle prop- 
erties. The spin resolved one-particle DOS calculated at 
T = for different values of U shows the expected insu- 
lating behavior with a clear gap at the Fermi energy for 
all U > 0. In particular for small values of U <C W, the 
DOS shows nicely developed nesting singularities at the 
gap edges, which qualitatively follow the predictions of 
a weak coupling theory>Ai With increasing U these fea- 
tures get more and more smeared out, and for U > W the 
spectra resemble those of the Mott insulator. 4 Note that 
the appearance of a gap in the DOS is of course accom- 
panied by a vanishing imaginary part of the one-particle 
self-energy in this region. Neither the development of the 
DOS nor the magnetization as a function of U shown in 
the inset to Fig. |3 provide any evidence as to whether 
the limits U <C W and U > W will be linked smoothly 
or by some kind of transition. 



B. Spin dynamics 



Since the ground state of our system is the Neel state, 
spin excitations require a minimum excitation energy, 
the spin gap A s , which conventionally is read off 
Smxi (Q, u)+iS) evaluated at the antiferromagnetic wave 
vector Q = (7r, 7t, ...). However, the gaps at other q 
vectors will be equal to or larger than A s . Thus, even 
after summing over all q-values, the size of the gap in 
3mxi( w + wm be determined by A s . The quantity 
$sTnx±{w + iS) on the other hand can easily be calculated 
from the NRG once the DMFT has converged. The re- 
sults for three typical values of U/W are shown in Fig. 01 
displaying a nice spin gap A s as u) — > 0. Evidently, the 
value of A s first increases with increasing U but then 
decreases again, as is to be expected from the mapping 
of the Hubbard model to an antiferromagnetic Heisen- 
berg model with J oc 1/U at large U. From the calcu- 
lated ^mx±(uj + iS) one can directly extract the values 
for A S (U). The results will be discussed below together 
with the charge gap obtained from the optical conduc- 
tivity (see Fig.[7J). 



Another interesting quantity is the dynamical mag- 
netic susceptibility, whose low-energy behavior gives fur- 
ther insight into possible differences in spin and charge 
dynamics. In principle, it is also possible to calcu- 
late this quantity as a function of wave vector q within 
the DMFTAiS However, this requires the calculation 
of the local irreducible particle-hole self-energy, 20 which 
is presently not possible within the NRG. Nevertheless, 
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FIG. 4: Imaginary part of the local transverse magnetic sus- 
ceptibility as a function of uj/W for U/W = 0.25, 0.5 and 
1, which shows a well-defined gap A a as u) — + 0. Only the 
part for uj > is shown. Note that A s first increases with 
increasing U, but eventually decreases again. 

for the current investigation, a reasonable approximation 
can be obtained from the local magnetic susceptibility, 



C. Optical conductivity and optical gap 

The optical conductivity resulting from the spectra in 
Fig. |21 is shown in Fig. [SI Apparently, the overall behav- 
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FIG. 5: Optical conductivity of the half- filled Hubbard model 
in the Neel phase for T = as a function of U. The full lines 
represent the calculated data, the dashed lines a fit with the 
function w-cr(w) = 9m |e l< ^ (u — luq + i r y)~ a } (see text). The 
inset shows the curves for U/W = 0.25 and U/W = 0.75 using 
a logarithmic scaling. 



ior seen in the DOS has its counter part in <j{u>). For 
small values of U, one finds a threshold behavior with 
a singularity, whereas for large U the optical conductiv- 
ity closely resembles the one found in the paramagnetic 
insulatoriiS Obviously, there are at least two interesting 
features in cr(cj). First, the behavior of cx(a>) in the vicin- 
ity of the maximum and second the actual value of the 
optical gap, i.e. the energy at which (j{lo) = 0. 
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In order to address the first point we adopt the follow- 
ing line of reasoning. In the Hartree limit, i.e. without an 
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FIG. 6: Dependence of the fit parameters ljo, a and 4> in @ 
on U . The lines are meant as guides to the eye. Note the 
rather well defined change in (a, cj>) from (a,cj>) = (l/2,7r/2) 
to (a, <f}) = (1, 0) around U/W = 0.75. 

imaginary part of the self-energy, an approximate evalu- 
ation of J3J yields 



ui ■ cr(u>) oc 



e(uj - 2A ) 
y/uj - 2A 



with Ao = Um s /2 and m s = (n-f — raj,). Since this be- 
havior is goverened by square root singularities in the 
integrand in (see e.g. the explicit formula derived in 
the appendix), it is reasonable to assume that for a finite 
imaginary part of the self-energy the above singularity 
will become an algebraic function 



u> ■ cr(u>) oc Sra 



,i<p 



(5) 



with a general exponent a. The quantity 7 approxi- 
mately cares for the finite imaginary part introduced by 
the one-particle self-energy and (f> allows for a more com- 
plex mixing of real and imaginary parts in the integral 
Q. The function (JSJ describes the behavior of <t(lu) in 
the vicinity of the maximum very nicely for all values of 
U (see dashed lines in Fig. EJ) ; note that from the inset to 
Fig. [5] it is evident that for small U this algebraic form 
has the tendency to overestimate the optical gap, while 
at large U it is clearly underestimated. 

Let us now turn to the behavior of the parameters 
luq, ol and shown in Fig. As U — > 0, we expect 
that ojq = 2Ao = Um s , a — 1/2 and <f> — tt/2, i.e. 
ui ■ a(u) oc 5Re (u> — ujo + iS) = Q(uj — ujq)/^/ui — ujq. 
We indeed find the anticipated square-root singularity; 
however, even for small U /W, the value of loq signifi- 
cantly deviates from the Hartree value, being systemati- 
cally smaller but obviously approaching it as U — ► 0. 



For values U > W, the behavior of u> ■ cr(ui) is best 
described by a Lorentian, which becomes apparent from 
the values of a and <f> obtained in this region, viz awl 
and 4> = 0, meaning u> ■ a(ui) cx Qrre (uj — ujq + 17) oc 
1/ ((w — ujq) 2 +7 2 ). In addition, the results for cjq to- 
gether with m s w 1 indicate that ujq w U, in agreement 
with the predictions of the Mott-Hubbard picture. 9 

The behavior of the optical gap A c , together with the 
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FIG. 7: Optical gap A c /2 (circles), spin gap A s (squares) and 
double occupancy (diamonds) as function of U. The inset 
shows the gaps scaled with Um a /2. Dotted lines are meant 
as guides to the eye. For small U both charge and spin gap 
are identical, while for large U we find A c oc U (see inset) and 
A s oc 1/(7 (full line in main panel). 

spin gap A s obtained from Qmx±{uj + iS) and the double 
occupancy (njnj), is displayed in Fig. The details of 
the method used to obtain A c are discussed in appendix 
B. 

For small U/W, the optical gap is exactly twice as large 
as the spin gap and, as becomes apparent from the inset, 
approaches the Hartree value m s ll as U — > 0. Again 
both quantities deviate systematically and by the same 
amount from the expected Hartree values even for the 
smallest U. Thus, even for U/W << 1 correlation effects 
are important and significantly modify the predictions 
from Hartree theorypiS^ For large U, on the other hand, 
we find A c oc U — W, consistent with Mott-Hubbard 
localized states; furthermore, A s oc 1/U as expected from 
the mean-field theory of the Heisenberg model with a 
Joe 1/U. 

We find, however, no evidence that the Slater limit at 
U/W -> and the Mott-Heisenberg limit at U/W -> 
00 are separated by some kind of phase transition. All 
results, including the variation of (n^n^) seen in Fig. 
rather indicate that a smooth crossover takes place for 
U/W w 3/4. 
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IV. CONCLUSION AND OUTLOOK 

While in the paramagnetic phase of the Hubbard 
model at half filling, when artificially extended to T = 0, 
a true phase transition from a correlated metal to a Mott- 
Hubbard insulator at U c ~ W has been established, the 
situation in the physically more relevant Neel state has 
not been investigated in similar detail up to now. As a 
first step into this direction the properties in the ground 
state of the Hubbard model at half filling with particle- 
hole symmetry have been discussed. We did confirm that 
the physical properties at small and large values of the 
Coulomb parameter U can be well described within a 
Slater and Mott-Heisenberg picture, respectively. In con- 
trast to the paramagnetic Mott-Hubbard MIT we could 
not find any solid evidence for a similar transition in the 
Neel state; our data rather suggest a smooth crossover, 
which occurs at a value U < W . Even the double occu- 
pancy, which in the case of the paramagnetic MIT is an 
indicator of a phase transition, does not show any sign 
of an anomaly here. 

A novel quantity we discussed was the spin gap, which 
we extracted from the local transverse spin susceptibil- 
ity. The general behavior and size agree very well with 
exactly known limits. This shows that at least in cases 
where a well defined spin gap exists that becomes mini- 
mal at special points in the Brillouin zone, even the in- 
spection of purely local dynamical susceptibilities can be 
sufficient. 

There are, however, still several unanswered questions. 
First, the analytic form of the optical conductivity close 
to the optical gap and the precise value of this gap could 
not be obtained at present due to numerical problems 
when evaluating the integral J3J. Especially for a more 
quantitative comparison with experiment this has to be 
improved in future work. Second, comparison with the 
data for V2O3 from ref. |5|show a nice agreement for the 
kinetic energy (obtained using the optical sum rule), but 
fail completely concerning the size of the optical gap. 
This is shown in Fig. [SI where we plot the kinetic energy 
scaled to its value for U — (circles and left scale) and 
&c/W (squares and right scale) vs. U/W. The open 
circle and square represent data for the kinetic energy 
and charge gap, respectively, extracted from ref. |^ for 
a V2O3 sample with Tjv « 50K. Obviously, with the 
present model, one largely overestimates the optical gap 
at intermediate values of U. 

Of course the present investigation did concentrate on 
the simplest situation, viz a system with perfect particle- 
hole symmetry. In reality, electron hopping beyond near- 
est neighbors will destroy antiferromagnetism at small 
values of U and consequently lead to different gaps at in- 
termediate values of U. On the other hand, the gaps at 
large U are controlled by Mott-Hubbard physics and will 
most likely change only little. A similar line of argument 
has in fact been invoked in ref. 0, too. A recent study 
of the properties of the magnetically frustrated Hubbard 
model indicates that this scenario is indeed very likely. 27 




U/W 



FIG. 8: Kinetic energy scaled to its value at U — (circles and 
left scale) and full optical gap A c /W (squares and right scale) 
vs. U/W. The open circle and square with errorbars represent 
values for the scaled kinetic energy and A c /W respectively, 
extracted from ref. |3 for the sample of V2O3 with Tjv ~ 50K. 
Note that in ref. 5 D = W/2 and A = A c /2. 



Evidently, a further detailed investigation of the opti- 
cal properties in the Neel state, in particular with more 
realistic band structures, is necessary and surely highly 
interesting. Work along this line is in progress. 
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APPENDIX A: FURTHER EVALUATION OF © 

In this appendix we present details of the further eval- 
uation of the energy integrals in equation Q for the hy- 
percubic lattice in the limit D — > 00. In that case the 
density of states becomes a Gaussian, i.e. we need to cal- 
culate the two integrals 



dee e 'A a (e,Lu')A s (e,uj'+uj) 



and 



/oo 
dee- t2 B a (e,oj , )B a (e,uj , + uj) . 
-00 

It is now convenient to split up the spectral functions 
into two parts, i.e. 

A a (e,uj) = A-(e,u) + A+(e,Lj) 



7 



with 
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and Co- as denned after equ. @. In the same way we 
write 



B (T {e,u) = B-{e,u)-B+{e,u) 



where now 



2tt VUC*~±t 



Using this notation and collecting equivalent terms, it 
can easily be verified that the following four integrals 
need to be evaluated: 

dee-<?B-{e,u')B-{e,u'+u) , 

-OO 

/OC 
dee-^B-^uj^B+^uj'+oj) , 
-00 



and 



h = 2 



h = 2 



dee'* A-{e,u)')Ar(e,u>'+u)) 



dee'* A-{e,J)A+(e^'+uj) 



The further evaluation will be demonstrated for the first 
term. Using the notation 



a = \/CtC;t and = ^ 'Q<tC& 

uj' -\-iS 



uj'-\-uj-\-iS 



we may write 

B-(e,Lj')B-(e,u'+u) = 
1 

~4^2 



a — e a — e 



0-e 0- 



where the bar above a term denotes complex conjugation. 
The terms inside the brackets can be expanded further 
to yield 
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a-0\a-e — eJ a — (3\a — e (3 — e 
1/1 1 \ 1/1 1 



a — (3\a — e (3 — e J a — 0\a — e (3 — e 
If we introduce the Faddeeva function 



i f°° e~* 
w(z) = — I dt 

7T J_m Z-t 



e z erfc (— iz) 



for complex arguments z with Smz > 0, we obtain 
h = 



1 

2tH 



w{a) — w{(3) w(a) — w{j3) 



a — (3 



a — (3 



w{a) + w{0) w(a) + w{(3) 



a — (3 



a- (3 



where we have made use of the relation u>(— z) = w(z). 
The remaining three contributions can be obtained in a 
similar fashion. Finally, combining complex conjugate 
expressions, we arrive at the following result for the four 
integrals: 



h = 



w(a) — w(0) w(a) + w{(3) 



a- (3 



a- (3 



T 1 - 

h = W71 



w{a) — w((3) w(a) + w{0) 



a + (3 



a + (3 



(Al) 



(A2) 



1 / w(a)-w(0) - w(a)+w(0) 

I 3 = - 3to 7(5 -yo = 

7T \ a - (3 a - 



1 j wja) - w(0) x w(a)+w(0) 

1a = — sm 70 = 70 

7T \' a + f3 a + 



where we have introduced 



(A3) 



(A4) 



and S = 

u'+iS vCcCff 



>'-\-uj-\-iS 



A further analytical evaluation of the remaining integra- 
tion over ui' in equation (QJ using eqs. (|A1|) - 1A4II is 
possible only for ( a — ► oj — ctAo + iS. In this case, the 
square-roots appearing in the functions a and lead to 
a typical threshold behavior of the form 22 



uj ■ <t(u>) (X 



e(to - 2A ) 
- 2A 



The appearence of this threshold singularity also shows 
that a further numerical evaluation of the remaining in- 
tegral over u' in equation Q will become problematic in 
regions where the imaginary part of the one-particle self- 
energy becomes small, because the integrand will develop 
a strongly singular behavior. In particular, this makes a 
precise numerical evaluation of the optical conductivity 
near the threshold impossible. 



APPENDIX B: THE OPTICAL GAP 

While the definition of the optical gap is straight for- 
ward, the extraction of numbers from the numerical data 
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appears to be rather problematic for two simple reasons. 
First, the spectra calculated with NRG have an unavoid- 
able intrinsic broadening, which becomes especially se- 
vere for the Hubbard bands at larger values of U/W. Sec- 
ond, as uj <C wo, the imaginary part of the one-particle 
self-energy becomes negligible, and the singular struc- 
ture of the integrand IjAlfl - l|A4(l entering into (@J makes 
additional broadening necessary to allow for a stable nu- 
merical integration. Together both effects very efficiently 
mask the true w-dependence close to the optical gap, in 
particular for larger U/W. 

In order to nevertheless have an unambiguous working 
procedure that allows to extract a reasonable approxi- 
mation to the true optical gap from our numerical data, 
we postulate that u> ■ <t(uj) cx Q(uj — 2Ao) • (u> — 2Ao) a for 
lu in the region where JSJ starts to deviate substantially 
from the data and choose a minimal a such that it pro- 
duces a reasonable fit for all values of U (see Fig. [3] for 
selected results) . We find a — 5 /22& and an optical gap 
A c which is consistent with the spin gap A s as U — > 0. 
The good agreement of these two differently calculated 
quantities (see inset to Fig. [7J also serves as an a poste- 
riori check for the procedure used to determine A c . In 
view of a possible comparison to experimental results 5 
this situation is, of course, not satisfying. For this pur- 
pose a more thorough and possibly analytical evaluation 



of a(uj) close to A c would be desirable. Unfortunately, 
the complicated form of the integrals in Q so far have 
allowed for an analytical evaluation only in the Hartree 
limit. 
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FIG. 9: Results for the fits to cr(ui) (full lines) with function 
© for uj w uj (dashed lines) and 9(w - 2A ) • (uj - 2A ) 5/2 
in the low-w region (dotted lines) for several values of U/W. 
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